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Summary. Probability models are proposed for passage time data collected in experiments with a de- 
vice designed to measure particle flow during aerial application of fertilizer. Maximum likelihood estima- 
tion of flow intensity is reviewed for the simple linear Boolean model, which arises with the assumption 
p '• that each particle requires the same known passage time. M-estimation is developed for a generaliza- 

tion of the model in which passage times behave as a random sample from a distribution with a known 
mean. The generalized model improves fit in these experiments. An estimator of total particle flow is 
constructed by conditioning on lengths of multi-particle clumps. 
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1. Introduction 



■ Measuring the outflow of granular particles from an airborne spreader during the aerial application 

£Nl ! of fertilizer or pesticide presents agricultural engineers with a difficult problem. The goal of uniform 

distribution over a targeted field requires knowledge about flow rate of the material as it is dropped 
from the aircraft. Wind s peed, air speed, granule properties, humidity, and temperature have been 
identified ( Casadv et al. (1993)) as factors which can lead to variability in these outflow rates and 



q 

hence amounts of material that reach the target. Typically, applicators are calibrated annually so 
that they achieve an average target flow rate. In practice, pilots use a simple lever-operated gate 
to change the flow rate in order to account for extreme values of these factors. This adjustment is 
based on intuition, without any feedback from measurement of particle flow. 



One approach to providing the pilot with more information uses an optical sensor device ([Grift and Hofetee 



■ ( 19971 )) which measures the velocity (in meters per second) and size of clumps of particles as they 



flow through the spreader duct. This device has two photo-sensitive arrays of optical sensors that 
receive a signal from a light source. As a particle passes an active area, it blocks this light thereby 
interrupting the signal received by the sensors. As long as all of the sensors in the array are receiving 
a high signal, the channel is classified as unoccupied and this is taken as an indication that there are 
no particles flowing through at that instant. If the signal to any one of sensors is interrupted, this 
is interpreted as the presence of at least one particle, constituting a clump, in flow. The two sensor 
arrays are 0.00078 meters apart and it is possible to measure the time in seconds that it takes a 
clump to move from one array to the other, Atf. The total time that either array is blocked, Atb, is 
also measured, facilitating calculation of velocity in meters per second, v — 0.00078/ Atf and clump 
length in meters, CL = vAtj,. These observable clump lengths, either in terms of physical length 
in meters or time in seconds, are the basis for inference ab out particle flow in the system. Such a 
me asurement device is cal led a type II counter (iPvkd (Il958h 1. 



Grift et al. ( 2001 ) and GriftJ ( 2002 ) carried out bench-scale experiments to evaluate the optical 
sensor device in situations designed to simulate the flow of fertilizer particles through an airborne 
spreader duct. In these experiments, a known number of spherical particles with a known mean 
diameter of 4.45mm was dropped from predetermined heights through a duct on which the sensor 
device was installed. The heights from which the particles were dropped was controlled at several 



2 



J.A. Osborne and T.E. Grift 



Run #20 



LO 

c\i 



o 
c\i 



Q 



LO 
O 



1 1 11 ITlWlTrTTHi^^— 



5 10 
Clump length(msec) 



15 



Fig. 1. Probability histogram of N(t) = 1790 clump lengths in msec 



values to simulate a range of particle velocities and flow rates. A histogram representing the dis- 
tribution of particle clump lengths in units of time, obtained from one run of these experiments, is 
shown in Figure 1. The relative frequencies for clump lengths (in msec) are based on dropping 4000 
spherical steel particles (actually BBs) from a fixed height. 

In this paper, simple linear Boolean models are used to describe the clump length data generated 
by the Grift experiments, thus providing a basis for inference about flow rates and total particle 
flow during a dispersal period. In particular, flow intensity is quantified by a single rate parameter 
in a simple Boolean model. Maximum-likelihood is reviewed in cases where particles require a fixed 
time for passage and an M-estimator is obtained in more general cases. Assessment of total particle 
flow utilizing this estimator is also developed. 

Section 2 introduces the Boolean model s and estab lishes notation and terminology. Results 
for the clump length distribution derived in Hall (19881) are used to develop M-estimation of flow 
intensity and the M-estimator is compared with maximum likelihood and other moment estimators 
by simulation. In section 3, two estimators of total particle flow are proposed, including one obtained 
by derivation of the conditional expectation of the number of particles in a clump given clump 
length under the equal diameters model. Simulations are carried out to give an assessment of the 
performance of this estimator in the random passage times model. The methods are evaluated based 
on their performance with the experimental data in Section 4. Section 5 concludes. 



2. Estimation 

To obtain a probability model for the clump-length data, particles are assumed to be identically 
spherical with a known diameter, do, and to arrive at the sensor according to a homogeneous Poisson 
process with unknown intensity A. Passage of particles is assumed to continue unabated upon arrival 
at the sensor. In one ver sion o f the model, the particles are travelling at a constant velocity, say vq, 
and the segment length ([Hall (1988)), or time required for any single particle to pass the sensor, is 
constant at to = do/vo. In a second version, velocities or diameters are assumed to vary in such a 
way that segment lengths behave as a random sample from a population with a known mean ji and 
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Fig. 2. A diagram of particle flow measurement using a type II counter 



an unknown variance a 2 . The two models will be referred to as deterministic segment length (DSL) 
or random segment length (RSL) models, respectively. 

Suppose that particle flow is observed for t time units. Let the number of particles arriving at 
the sensor in this time period be denoted by A(t). Let N(t) denote the number of complete particle 
clumps observed by time t. Let Yi,Yi, ■ ■ ■ , ^Ar(t) denote the lengths of these clumps and Z±, Z2, ■ ■ ■ 
the spacings between them. Let the unobservable number of particles comprising clump i be called 
the clump order and be denoted by iQ. 

Figure 2 illustrates the clumping process using an example with A(t) = 7 particles arriving at a 
sensor at times 1.9, 5.9, 6.8, 7.5, 11.6, 12.8 and 17.1 msec during an observation period of t = 20 msec. 
If particles are assumed to have diameter 4.45 mm and to be travelling at a constant velocity of 
2.225 mm / 'msec, the passage time required for each, or deterministic segment length, is do = 2 msec, 
leading to four clumps of lengths y\ = 2, 1/2 — 3.6, j/3 = 3.2, j/4 = 2 msec that exit the sensor at times 
3.9, 9.5, 14.8 and 19.1 msec, respectively. Spacings between clumps would be of length z\ = 1.9, Z2 = 
2.0, Z3 = 2.1 and Z4 = 2.3 msec and the four clump orders would be k\ = 1, = 3, ks = 2, kn = 1. 

The particle clumps constitute a coverage process on one dimension. lHall (|l988l) describes 
the process as a simple linear Boolean model; simple because the clump-lengths are line segments 
and linear because the events occur in one dimension, the time line. Linear Boolean models also 
arise as linear transects from higher dimensional convex-grain Boolean models. In the language of 
queueing theory, the number of particles in a clump at the sensor at a given time forms an M/D/oo 
queue in the DSL model and an M/G/co queue in the RSL model and clum p-lengths are called 
busy periods. There is much literature on these models from queueing theory ( Daley ( 200lh ) . For 
statistical inference for the distribution of diameters or more complex quan ti ties descr i bing t he grain 
process, or for Boolean models in higher dimensions, see iMolchanovl (jl997h . lHandievI (|l999f l derived 
a discrete approximation to the distribution of clump-length in the linear Boolean model and used it 
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for likelihood inference. ICrespi et al. I (|2005l) have employed the linear Boolean model for monitoring 
events of viral activity in human subjects. 



2.1. Likelihood 

Specification of the clump-length density, f(y\ A ) is di fficult outside of the case where particle di- 
ameters are degenerate. In the DSL model, iHal (|l988f ) has shown that the density has point mass 
to, and is otherwise given by 



e~ At ° at y 



/(y;A,to) = A- 



-Ai 



E 



(-i)j 



{X(y - (j + l)t )p- 1 e^ At HA(2/ - (J + l)t„) + j} 



where y > to and s is the largest integer such that to < y/s. The continuous part of the density is 
uniform over (to-, 2io), and decreasing for y > 2to- For small A, f(y; A, to) can be approximated by the 
uniform distribution on (to, 2to), for large A it can be approximated by an exponential distribution. 
The fitted density f(y; A = 0.40, to = 2.00msec) overlays the probability histogram of experimental 
clump-lengths in Figure 1. 

For RSL models, likelihoo d inference is difficult because of the complexity of the clump-length 
distributions (|Handlevl (|2004h ). In the DSL model, an approximate likelihood function can be 
specified by ignoring the residual lifetime of the process. The residual lifetime is the duration of the 
last incomplete clump or spacing. A clump is a singleton if there are no arrivals within to time units 
of the start of the clump, an event which occurs with probability e _At °. Let M\ = : yi = to} 

denote the number of singleton clumps. By independence of clump-lengths, the approximate partial 
Boolean likelihood can be factored into components for singleton point masses and multi-particle 
clump-length densities: 



£(A;yi, . . .,y N (t)) 



-mi A^o 



singletons ' :w>t ° 



multi-particle lengths. 



Spacings z\, 22, • • • are not available for the experiments analyzed in section 4. For cases where the 
Zi are available, an approximate complete Boolean likelihood may be obtained by multiplying the 
partial likelihood by the likelihood from an exponential random sample, \ N e~ x ^ Zi . 

For large t, the maximum likelihood estimator of A based on £ is approximately normally dis- 
tributed. However, the analytic expression for the Fisher information is unwieldy, particularly for 
large clump- lengths, where the degree of the polynomial components of the clump- length density 
is high. Alternatively, approximate confidence regions can be constructed from the likelihood ratio 
test statistic, which has an approximate \ 2 distribution on 1 degree of freedom. 

In the RSL model, where segment lengths are distributed as a random sample from a known 
distribution with distribution function G(x;9), the clump-length density and resulting likelihood 
are considerably more complex. Let fnsL(y, A, 9) denote the clump-length density, which depends 
on the unknown parameters, A and 0. Ignoring the residual lifetime, the partial likelihood of the 
complete clumps is then 
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Halll (|l988l ) shows that the Laplace transform 7 of Y is 



7(s) 



- A 



, / cxp{-st-A / {1- G(x;0)}dx}dt 
Jo Jo 



Stadjel ( 1985f ) obtains the clump- length distribution function Frsl{ii) by inversion of 7, but it is an 



infinite sum of self-convolutions of a function that may involve an integral with no analytic solution, 
making inference based on Crsl difficult. 



2.2. M-estimation 

An important issue in estimation of A is robustness under model misspecification. Inspection of 
the clump-lengths from the experimental data, such as the run depicted in Figure 1, reveals that 
the number of clumps with lengths slightly in excess of to is greater than expected, so that the 
distribution between to and 2to is not uniform. This can be caused by variability in diameter or 
velocity or by errors of measurement. A desirable property for any estimator is robustness to this 
departure from model assumptions. 

For mean segment length fj, 7 the mean clump-length is given by 

eV _ 1 



in either the DSL or RSL model, regardless of the distribution of segment lengths (Hall (1988)). For 
known /x, consider the M— estimator A which satisfies 



e^ 1 - 1 

y = 



A 

A solution exists by the mean value theorem with E(Y; A) increasing in A. Though there is no 
analytic solution, the equation can be solved rapidly usi ng any root-finding pr o cedur e, such as the 
uniroot function in the R statistical software package ()lhaka and Gentleman (|!996t) L A starting 



point that works in simulations is given by A = (y — /x)/(2/i 2 ), which is the solution obtained using 
a second order expansion of e A/i about 0. An interesting aspect of the sampling distribution of A 
is that it is negative whenever y < /i, an event whose probability is small as long as Xfi is not too 
small. 

This estimating equation for A can be written 



5>(y*,A) = 



where ijjjy. A) = y— A 1 (e A<5 — 1). Large-sample theory for M-estimators, (see, e.g. Stefanski and Bo"osl 



can be used for inference about A. For a random sample of n clump-lengths y\, . . . , y n , the 
asymptotic distribution of A is given by 

Vi(A-A) N(0,C/B 2 ) 
where B and C are functions of A defined by 

B(X) = Ei-^iY^X)) 
C(X) = E^iY^X)). 
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Since ip is linear in Y, the expectation operations are straightforward: 

e^iXfi- 1) + 1 



B(X) 



A 2 



C(A) = Var(F;A). 

The variance of Y depends on the distribution of segment lengths. In the DSL model with to = M; 

Var(F) = \- 2 {e 2Xt * - 2X^e Xt * - 1). 

In the RSL model with segment lengths distributed according to the general distribution function 
G(x), clump- lengths have variance 



Var(F) = 2A-V M j (^exp X 



(1 - G{x))dx 



-\\dt-\-\e 



I) 2 



which can be estimated using the sample variance of clump- lengths, s y . Estimators for the variance 
of A are then given by 

_ 1 A 2 (e 2 ^ - 2Ve^ - 1) 



Var(A) = n 



(e^(A/i-l) + l) 2 



in the DSL model and 



Var Q (A) = n- 1 - 



{e**{\p - 1) + l) 2 

in either the DSL or RSL model. In large samples, approximate confidence intervals for A can be 
constructed from these estimates along with the normal approximation for A. 



2.3. Other estimators 

For the DSL model with common deterministic passage time to, other method-of- moments (MOM) 
estimators can be constructed using only the clumpcount (N(t)) and singleton count ( Mi) statistic s. 
The sequence of i.i.d. sums {Zi + Yi} is a renewal process. Elementary renewal theory ( Coxl ( 19621 )) 
yields that as t — > oo, 

N(t) - tlu R c 

V } N(0, 1) 

where and a\ denote the mean and variance of a randomly sampled renewal period. In DSL 
model with deterministic common passage time ioj 



(j, R = E(Z + Y) = A _1 e At ° 

CT 2 = Var(Z + F) = A~ 2 (e 2A *° - 2At e At °). 



Moments for N(t) are then 



-At 



E[N(t)] « Xte" 
Var[JV(i)] w Xt (e" At ° - 2At e _2At °) . 

The probability that a randomly selected clump is a singleton is e _At ° so that E(M±) = Ate~ 2At ° 
A MOM estimator based on the observed number of singletons is then 
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Grift et all (1200 lh and iGriftl (|2002h base estimation of total mass now on this estimator. Other 
estimators of A can be co nstruc t ed by consideration of vacancy, V ~ ^ Zi, or total time that that 
the sensor is unoccupied. Halll (|l988l ) develops asymptotic theory for a number of vacancy-based 
estimators. Measurements of V were not available from the experiments discussed in section 4, and 
vacancy-based estimators are not considered further. 



2.4. Simulation 

Simulations were undertaken to provide some information about the performance of these estimators, 
with three goals in particular: a comparison of the efficiency of the moment estimator A relative to 
the MLE under the DSL model, an investigation of the robustness of the MLE under the RSL model 
and a comparison of coverage probabilities of confidence intervals resulting from the two variance 
estimates of the asymptotically normal M-estimator, A. Particle arrivals were generated according 
to a Poisson process. Three cases with an increasing degree of clumping were simulated using flow 
intensities of A = 0.1, 0.2 and 0.3. Two times were considered for the length of the total observation 
period, t — 1000 and t = 10000. Preliminary experiments with particles far enough apart so that 
there was no clumping indicated that measured passage times were normally distributed. So, passage 
times for individual particles were generated from a normal distribution with a mean of /i = 5 with 
three different standard deviations, a = 0, 0.5, 1. The first of these standard deviations leads to the 
DSL model, the others to RSL models. The approximate mean clump counts for the DSL model 
were E[K] « 1.6,2.7,4.5 for the three flow rates, A = 0.1,0.2,0.3, respectively. The simulation 
experiment then had a crossed 3x2x3 design, with n = 500 independent datasets generated per 
combination of A, t and a. Normal plots and Kolmogorov-Smirnov statistics did not indicate any 
obvious non- normality for either the MLE or A. 

Table 1 summarizes the results of the simulation. The bias of the M-estimate relative to A and 
the efficiency relative to the MLE are given in the middle section. Though the bias of the MLE 
formulated under the DSL model dissipates with increasing A or t, it does not exhibit robustness 
to heterogenous segment lengths, in the sense that it has larger variance than the M-estimate. 
Empirical coverage probabilities for 95% confidence intervals based on the LRT and those of the 
form A ± 1.96SE where SE denotes the appropriate estimated asymptotic standard error from 
Section 2 are given in the right section of Table 1. For the shorter simulations (t = 1000), there 
is a tendency for coverage probabilities based on A to be low. For datasets with a larger number 
of clumps (t = 10000), the nominal coverages for intervals based on A are reached. With n = 500 
simulations, the Monte Carlo standard error is such that any sample proportion less than 0.934 
is significantly less than the nominal 0.95 with comparisonwise error rate 0.05. Additionally, the 
intervals around the M-estimate that use the standard error, SEg, which is a function of the sample 
variance of the clump-lengths, appear to do better for the RSL models with large N(t), particularly 
for the noisy segment length a — 1 case. The likelihood ratio interval gives coverages consistent with 
nominal levels in simulations with the DSL model, but breaks down under the RSL model where the 
likelihood is misspecified. In summary, the recommendation based on these simulations is that the 
M-estimator is reasonably efficient under the DSL model and robust to the conditions of the RSL 
model. Confidence intervals based on the standard error SEq meet nominal coverage probabilities 
in large samples under either model. 



3. Estimation of total particle flow 

In the case where either A is known or variance in its estimation is negligible, total particle flow 
may be estimated by i?[A(t)] = At. When t is not available, another estimator can be formed by 
substitution of t w X^^i + 2 E(Zi) into the expression giving i?L4(<)] = A^ Yi + N(i). 
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Table 1. Simulation: relative efficiency and coverage probability of A 
estimators 



Parameters Rel. Rel. Coverage probabilities 



O" 


t 


A 


N(t) 


Bias 


Eff. 


LRT 


SE{\) 


SE G {\) 





1000 


0.1 


60.1 


0.01 


0.89 


0.966 


0.950 


0.944 





1000 


0.2 


73.2 


-0.01 


0.95 


0.956 


0.956 


0.934 





1000 


0.3 


66.7 


-0.01 


0.98 


0.948 


0.946 


0.944 





10000 


0.1 


605.8 


0.00 


0.83 


0.940 


0.942 


0.942 





10000 


0.2 


734.5 


0.00 


0.88 


0.960 


0.956 


0.950 





10000 


0.3 


669.6 


0.00 


0.97 


0.942 


0.938 


0.942 


0.5 


1000 


0.1 


60.2 


0.56 


8.3 


0.128 


0.922 


0.916 


0.5 


1000 


0.2 


73.0 


0.15 


2.2 


0.768 


0.910 


0.902 


0.5 


1000 


0.3 


66.3 


0.05 


1.1 


0.940 


0.952 


0.950 


0.5 


10000 


0.1 


606.1 


0.57 


82.4 


0.000 


0.942 


0.946 


0.5 


10000 


0.2 


736.1 


0.15 


17.2 


0.002 


0.946 


0.954 


0.5 


10000 


0.3 


670.1 


0.05 


3.3 


0.646 


0.938 


0.940 


1 


1000 


0.1 


60.5 


0.56 


7.3 


0.136 


0.914 


0.940 


1 


1000 


0.2 


73.5 


0.13 


1.9 


0.798 


0.922 


0.920 


1 


1000 


0.3 


66.9 


0.03 


0.98 


0.952 


0.944 


0.942 


1 


10000 


0.1 


605.3 


0.56 


71.9 


0.000 


0.922 


0.952 


1 


10000 


0.2 


735.7 


0.14 


13.6 


0.016 


0.924 


0.946 


1 


10000 


0.3 


668.6 


0.04 


3.0 


0.694 


0.944 


0.954 



In the DSL model, clump orders (Ki,K2, ■ ■ .) may be shown (IPippenger to be geomet- 



rically distributed. A clump is of order one (Ki = 1) if there are no arrivals within to time units 
of the start of the clump, which occurs with probability e~ At °. A clump is of order two if there 
is exactly 1 arrival within to units and none in the next to time units, an event which occurs with 
probability (1 — e -^t ^ e ~\t an( j SQ on _ _ _ are ^ independent geometric random variables 

with support on positive integers: 

Pr(Ki = k) = (1 - e -^k-i e -\t for = 1, 2, . . . 

with E(Ki) — e xt ° and Var(iQ) = e 2At ° —e xt °. If the system is vacant when observation ends at time 
t, then total particle flow may be expressed as the sum of these clump orders: A(t) = K\+- ■ -+K N ( t y 
If the system is occupied at time t, there is a partial clump that contributes a relatively small amount 
of particle flow for large t. Expressing total particle flow A(t) as the sum of clump orders each 
with mean e xt ° suggests the estimator Ai(t;X) = N(t)e xt °. When evaluated at the M-estimator 
A, with mean passage time /i = to, the two estimators of total particle flow become equivalent: 
Ao(X) = N{t)e x » = A 1 {t;~X). 

More efficiency might be gained by conditioning on the clump lengths. The estimator p(y) of an 
individual clump order which is a function of the clump length y and minimizes the mean squared 
error E[(K — p(y)) 2 ], is the Bayes estimate, or p(y) — E(K\Y = y). An estimate of mean total 
particle flow E[A(t)] — E\^2 Ki] is then given by summing over clumps: 

N(t) 

A B {t;X) = Y J E{K l \Y l ;\). 

i=i 



Of course E[Ki\Yi = to; A] = 1. The approach used by Haj (l988) to derive the clump length density 



f(y) in the DSL model may be extended to obtain the conditional mean of clump orders, E(K\Y). 
Let the beginning of a clump be the origin and let k denote an integer greater than unity. The joint 
event K — k and Y € (y, y + dy) occurs if and only if there is a particle arrival at {y — to, y — to + Ay), 
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no arrival in (y — to + Ay, y), exactly k — 2 arrivals in (0, y — to), and the nearest neighbor of each 
of these k — 2 arrival times is not further than to time units away. Since the first three of these 
conditions are independent and the fourth is conditionally independent of the first two given the 
third, the joint probability of these four events is the product 

\Aye- Xt ° {X{y ~ t0 Jt~\ - X ^p k -2 ( 



(fc-2)! "-'\y-to 

where p n (u) denotes the chance that the largest division formed by a random sample of n points 
taken from the unit interval does not exceed u. This probability is given by 



[u ~ 1] / +1 \ 

Pn (u) = E " 7 ) {1 - ju)n 



3=0 

= l-(„ + l)(l- u )»+( n + l ) (l-2«)"-... 

where [•] denotes the largest integer not exceeding the argument. Division by f(y) and differentiation 
with respect to y yields the conditional density 

Pr(K = k\Y = y) = T . --. p k -2 ' 



f(y) (fc-2)! ^-'\y-to 

If s — [y/to], then summation over positive integers yields an exact expression for the conditional 
mean: 

oo 

E(K\Y = y) = E kPr(K = k\Y = y) 

k=s+l 

Ae- At " ^ J\{y-to)) k - 2 to 

- Ae " t0 V fc ^-^))" 2 Vr iy( k - l \(i * 

Inspection of Pr(K = k\Y = y; A) reveals that for t < y < 2t , K has the translated Poisson 
distribution with mean and variance that are linear in y. For larger y, numerical evaluation of 
E(K \Y = y) be difficult. Inspection of plots for larger y and various values of A indicates that after 
a jump discontinuity of Atoe" At °(l — e _At °) _1 at y = 2t , approximate linearity extends to y > 2t - 
For cases where N(t) is large and there is heavy clumping, E(K\Y = y) can be approximated by 
linear interpolation to save computational effort. 



3.1. Simulation 

The performances of these estimators of mean total particle flow are compared using the simulated 
data from section 2. Error for either Ai or is, as a percentage of the mean particle flow is assessed 
using the relative root mean squared error, RRMSE: 



RRMSE(A(t)) - =1= 500-iJ2(Mt) - Mt)) 2 

A \t) y i 

where i indexes the 500 simulated datasets. Table 2 summarizes relative bias and RRMSE of 
estimates obtained by substitution of the M— estimates A into the expressions A\(t) = N(t)e XfJ - t and 
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Table 2. Error in estimation of total particle flow from 
simulations. 

Parameters Relative bias RRMSE 



(J-t 


t 


A 


Ai 


A B 


Ax 


A B 


o 


1000 


0.1 


-0.008 


-0.007 


0.042 


0.033 


o 


1000 


0.2 


-0.012 


-0.011 


0.048 


0.045 


o 


1000 


0.3 


-0.014 


-0.013 


0.049 


0.048 


o 


10000 


0.1 


-0.001 


-0.001 


0.014 


0.01 





10000 


0.2 


0.000 


0.000 


0.015 


0.014 





10000 


0.3 


-0.002 


-0.002 


0.016 


0.015 


0.5 


1000 


0.1 


-0.003 


0.178 


0.048 


0.185 


0.5 


1000 


0.2 


-0.008 


0.061 


0.051 


0.076 


0.5 


1000 


0.3 


-0.011 


0.013 


0.054 


0.052 


0.5 


10000 


0.1 


0.000 


0.184 


0.014 


0.184 


0.5 


10000 


0.2 


-0.001 


0.067 


0.015 


0.068 


0.5 


10000 


0.3 


-0.003 


0.022 


0.015 


0.026 


1 


1000 


0.1 


-0.002 


0.189 


0.058 


0.198 


1 


1000 


0.2 


-0.010 


0.063 


0.053 


0.079 


1 


1000 


0.3 


-0.011 


0.018 


0.052 


0.052 


1 


10000 


0.1 


0.001 


0.192 


0.018 


0.193 


1 


10000 


0.2 


0.000 


0.072 


0.017 


0.074 


1 


10000 


0.3 


-0.001 


0.026 


0.016 


0.031 



Asit'jX) for each simulated experimental condition. The estimation based on clumpwise estimated 
clump orders Ab, is competitive under the DSL model (cr t = 0) for smaller sample sizes, (t = 1000). 
It suffers from some positive bias in RSL models that appears to decrease as flow rate A increases, 
though it remains inferior to Ai despite smaller variance and higher correlation with A(t). In the 
RSL model, many singleton clumps have clump lengths slightly in excess of the mean singleton 
passage time \it and so have estimated orders in excess of 1. This may lead to a positive bias for the 
clumpwise estimators which is particularly acute when support is high near Y — fi t - This theory is 
supported by the poor performance under light clumping, when A = 0.1 and density near Y = [Xt is 
highest among values of A considered in the simulation. 

In summary, for minimal relative error, these simulations suggest the use of the simple A\{t) 
estimator, which is unbiased and involves less computation than the clumpwise estimator ^4s(t). A 
slight loss of efficiency under the DSL model may be offset by the superior performance in the RSL 
model. Expressed relative to total particle flow, the root MSE was not larger than 5.8% in any of 
the conditions simulated here. 



4. Experimental data 



An op ti cal se n sor w as used to measure clump lengths and clump velocities in experiments ([Grift et al 



d200lh : iGriftJ (j2002h ) m which a known number of spherical particles was dropped through a device 



simulating an aerial spreader duct. Various quantities of several kinds of particles (BBs, urea 
fertilizer) were dropped at several velocities. The data considered here include 10 runs with 4000 
identical steel particles (BBs) dropped from each of two heights and 5 runs with 2000 BBs dropped 
from a fixed height. Mean (y) and variance (Sy) of physical lengths (in mm) appear in Table 3 along 
with other statistics from the experiments. Division by mean velocity (v ~ 2.23mm/msec) was used 
to transform the measurements to the time line (in msec) to obtain Figure 1. In general, velocity 
was reasonably constant within a run of the experiment. 

The data were imperfect and some outlier removal was undertaken. For example, the counter 
returned several clumps with negative velocities or negative physical lengths, or sometimes both. 
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Table 3. Estimation from experiments with BBs 



Run 


N 


y 


Sy 


\[SE) 




A n 

^-B 


\ 


90^8 


5.22 


O.UI 


0.070 


(0.003) 


4041 


4921 


2 


9cnn 


5.22 


9 Q1 


U.U/U 


(U.UU3J 




4946 


•> 
o 


98Q1 


^ 9fi 


o.oy 


0.073 


(0.003) 


zinns 


4874 


4 




5.22 


o.uu 


0.070 


(0.003) 


dnnn 


4944 


K 
O 


9QQ0 


^ 1 fi 


9 88 


U.Uoo 


(U.UU3J 


"^QSfi 


4941 


fi 
U 


2941 


^ 90 


^ nn 


U.Uoo 


(U.UU3J 




488^ 


7 




^ 1 'l 


2.84 


n nfi/i 

U.U04 


(^U.UUo ) 


^QfiQ 


4Qnn 

iyuu 


8 


2956 


5.16 


2.90 


U.UOO 


(U.UUo ) 


3952 


4846 


9 


2894 


5.24 


3.12 


U.U / 1 


(U.UUo J 


3976 


4831 


10 


2914 


5.25 


3.11 


U.U / O 


(U.UUo ) 


4025 


4931 




1821 


fi 7fi 


1 1 77 




(U.UUo ) 


oyoo 


49QQ 
iiyy 


12 


1 77D 


fi 8^ 


1 1 ^fi 


n 1 89 




oy i u 


4^^ 


13 


1805 


6.80 


12.30 


0.179 


(0.005) 


4000 


4321 


14 


1748 


6.96 


12.57 


0.188 


(0.005) 


4038 


4333 


15 


1800 


6.85 


12.13 


0.182 


(0.005) 


4040 


4340 


16 


1784 


6.93 


14.82 


0.186 


(0.005) 


4089 


4403 


17 


1772 


6.93 


12.56 


0.187 


(0.005) 


4064 


4341 


18 


1788 


6.89 


13.28 


0.184 


(0.005) 


4052 


4346 


19 


1812 


6.78 


12.01 


0.178 


(0.005) 


3995 


4317 


20 


1790 


6.84 


11.98 


0.181 


(0.005) 


4005 


4330 


21 


746 


7.54 


17.10 


0.219 


(0.008) 


1981 


2143 


22 


791 


7.24 


13.20 


0.204 


(0.007) 


1959 


2141 


23 


777 


7.46 


13.15 


0.215 


(0.007) 


2024 


2184 


24 


774 


7.30 


13.23 


0.207 


(0.007) 


1941 


2102 


25 


745 


7.57 


13.39 


0.221 


(0.007) 


1989 


2133 



Additionally, each run contained a very small number of extremely short clumps, much less than 
the particle diameter, possibly due to matter other than the particles of interest blocking the sensor. 
The number of questionable clump measurements that were removed did not exceed 1% for any of 
the 25 runs. 

In these experiments, total particle flow is fixed and total flow time varies with run and is not 
observed. The opposite is true for the application of mass flow measurement during aerial application 
of fertilizer particles. The theoretical results regarding inference for the random particle flow A{t) 
for fixed t do not necessarily hold under the conditions of the experiment, where A(t) is fixed and 
t varies and is not observed. However, Table 3 provides some indication that estimates for total 
particle flow, A(t), have good empirical performance when it is treated as random, at least in these 
experiments. 

The observed value of the estimator Ai is given in the penultimate column. It appears to perform 
reasonably well under these conditions. The average of A\ over runs 1-20 is 9 and the root mean 
squared error from 4000 is 35.5, which is 0.9% of the target. There is some evidence of positive bias 
in the high intensity runs 11-20. A two-sided i-test of the hypothesis that -E[^Li] = 4000 under the 
conditions of runs 11-20 yielded a p-value of 0.065 on df = 9. 

Higher flow rates lead to more clumps per particle, fewer singletons, and larger variance in 
estimation of clump order, cither conditionally as in Ab or unconditionally, as in A\. The standard 
deviations of A\ under the light (runs 1-10) and heavy (runs 11-20) clumping conditions with 4000 
BBs were sj = 26.4 and Sh = 37.1, respectively. The estimates Ab, which are based upon the DSL 
model, exhibit substantial positive bias, as they did for data simulated under the RSL model. The 
same is true for the MLE of A. 

To assess the goodness of fit of the linear Boolean models, probability histograms of the clump 
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length data were checked for agreement with the estimated density f(y;X,d). One such check 
appears in Figure 1, which exhibits reasonable fit except for slightly lowered mass at the mean 
segment length, \i = 2msec and slightly more observations just above the mean segment length than 
expected under uniformity of this part of the density. All of the other histograms exhibited the same 
three distinctive features of a spike near this fixed segment length, near uniformity between one and 
two of these lengths and a long right tail. Quantile plots and Kolmogorov-Smirnov goodness-of-fit 
tests, for estimates in runs 1-10 or runs 11-20 do not indicate any non-normality in the distribution 
of AT, A or ii. 

5. Conclusion 

Two versions of a simple linear Boolean model are proposed to describe passage times of clumps 
of particles in a type II counter system; one assumes deterministically equal passage times for 
all particles, while the other assumes these to be distributed about a known mean with unknown 
variance. An M-estimator of flow intensity is developed that is intuitively sensible, computationally 
feasible, and robust to conditions where either particle velocity and/or diameters have substantial 
variability or are being measured with error by the type II counter. 

For total mass flow, A(t), two estimators are developed. The first is simply product of the number 
of clumps, N(t) and the estimate of the mean number of particles per clump. The second more 
complex estimator is the clumpwise sum of conditional mean clump orders (K), given clumplcngths 
(Y). In models where segment lengths were deterministically equal, the Bayes estimator exploiting 
the conditional mean clump order had relative root mean squared error not exceeding 5.0%, and 
always lower than that of the simpler estimator based only on the M-estimate of flow rate and the 
number of clumps. Under the most favorable conditions, with light particle flow and a long dispersal 
period, the relative error was as small as 1%. While the Bayes estimator did well in data simulated 
from the DSL model, it was outperformed by the simpler estimator in simulations where the segment 
lengths vary according to a normal distribution and in the bench-scale experiments. The relative 
root mean square error when using the estimator of based ranged between 1.4% and 5.1%. The 
relative root mean square errors for the experimental data were 0.6% and 1.1%, in the low and high 
intensity runs with 4000 BBs, respectively and 1.8% in the runs with 2000 BBs. 
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